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Observation of two distinct d xz /d yz band splittings in FeSe 
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We report the temperature evolution of the detailed electronic band structure in FeSe single¬ 
crystals measured by angle-resolved photoemission spectroscopy (ARPES), including the degeneracy 
removal of the d xz and d yz orbitals at the Y/Z and M points, and the orbital-selective hybridization 
between the d xy and d xz / yz orbitals. The temperature dependences of the splittings at the T/Z 
and M points are different, indicating that they are controlled by different order parameters. The 
splitting at the M point is closely related to the structural transition and is attributed to orbital 
ordering defined on Fe-Fe bonds with a d- wave form in the reciprocal space that breaks the rotational 
symmetry. In contrast, the band splitting at the T points remains at temperature far above the 
structural transition. Although the origin of this latter splitting remains unclear, our experimental 
results exclude the previously proposed ferro-orbital ordering scenario. 

PACS numbers: 74.70.Xa, 74.25.Jb 


Several experimental studies report the breakdown of 
the rotational symmetry in parent and underdoped com¬ 
pounds of Fe-based superconductors (FeSCs) [I}@] that 
is commonly refereed to as nematicity. Its origin is highly 
debated since both magnetic (MU and orbital [9in~2l fluc¬ 
tuations or orderings can lead to nematicity. Although 
strong support is given to magnetic-driven nematicity in 
iron-pncitides m where the orthorhombic lattice distor¬ 
tion is always accompanied by a collinear magnetic order 
at a temperature equal to or below the lattice transi¬ 
tion temperature, this mechanism is questioned in FeSe, 
which exhibits an orthorhombic lattice distortion below 
the distortion transition temperature T s ~ 90 K without 
any trace of magnetic order. As a direct signature of the 
electronic anisotropy between the x and y directions in 
the nematic state, previous angle-resolved photoemission 
spectroscopy (ARPES) studies [T4HT6] revealed a split¬ 
ting between the otherwise degenerate Fe 3 d xz and Fe 
3 d yz orbitals at the M point of the Brillouin zone (BZ). 
This splitting is widely believed to be a key evidence for 
ferro-orbital ordering in the nematic phase HUMS. 

In this letter we report the existence of two distinct 
splittings between the d xz /d yz bands of FeSe single¬ 
crystals near the Fermi level (Ep). We show that a first 
splitting decreases with temperature increasing, and dis¬ 
appears at about 100 - 120 K, which is slightly higher 
than T s , suggesting that it is caused by short-range or¬ 
bital order or fluctuations related to the structural transi¬ 
tion. This splitting has a d -wave form breaking rotational 
symmetry [20] that is the largest at the M point and 
that is inconsistent with ferro-orbital ordering. In addi¬ 
tion, we observe a splitting at the Y point that is rather 


insensitive to temperature up to 150 K, way above T s . 
Due to the strong orbital-selectivity of the hybridization 
between the d xy band and the d xz and d yz orbitals, we 
conclude that the splitting at the Y point is not simply 
due to spin-orbit coupling (SOC). 

High-quality single-crystals of (3 -FeSe were grown by 
the KCI/AICI3 chemical vapor transport method [2T] . 
The T c was determined to be 9 K from magnetization 
measurements and a structural transition is observed 
around 90 K. ARPES measurements were performed at 
the Dreamline beamline of the Shanghai Synchrotron Ra¬ 
diation Facility (SSRF) using a VG-Scienta D80 elec¬ 
tron analyzer, and at the Institute of Physics, Chinese 
Academy of Sciences, using a R4000 analyzer and a he¬ 
lium discharge lamp. The angular resolution was set to 
0.2°. Clean surfaces for the ARPES measurements were 
obtained by cleaving the samples in situ in a working 
vacuum better than 5 x 10“ 11 Torr. In the text, we label 
the momentum values with respect to the 1 Fe/unit cell 
BZ. 

We show in Fig. [l] the electronic band structure of 
FeSe, below the structural transition. The Fermi surface 
(FS) (Fig. 0 a)) is formed by one hole pocket centered 
at Z and two electron pockets centered at A. Based on 
local density approximation (LDA) calculations (Supp. 
Part I), we attribute the two elliptical electron pockets 
at A to d xz / yz bands in different twin domains, while the 
d xy electron pocket is not observed. A schematic repre¬ 
sentation of the FSs (T > T s ) and their areas are shown 
in Supp. Part III. As shown in Fig. [lj (5), the k z dis¬ 
persion along Y-Z is non-negligible, in agreement with a 
previous report [22] . The Fermi wave vector (kp) near 
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FIG. 1. (a) FS mapping along the Z(0, 0, 7r)-A(0, 7r, 7r) 

direction, recorded with unpolarized He la: photons. (6) FS 
mapping along the Y (0, 0, 0) - Z direction obtained in the 
a geometry. The red solid lines indicate the k z dispersion, 
(c - e) Band structure at T, Z and M (0, 7r, 0) along T-M 
or Z-A with T < T s , recorded in the cr geometry. (/ - /i) 2D 
curvature of (c - e). The red, blue and black lines indicate the 
dispersions of the d yz , d xz and d xy orbitals, respectively. The 
thick and thin lines in ( h ) correspond to 2 different domains. 
(i - k) EDC plots of (c - e). The black dots indicate the EDC 
peaks. 


the T point is ^ 0 . 077 r/a, while it is ^ 0 . 147 r/a at the 
Z point. The small kp can be clearly resolved from the 
cut at T displayed in Figs. 0c), 0/) and [ljz). Besides 
the d yz band, we also resolve a steep d xz and a flat d xy 
bands below Ep . Interestingly, the top of the d xz and d yz 
bands do to coincide, in contrast to LDA calculations but 
in agreement with a previous ARPES report m • From 
the energy distribution curves (EDCs) and curvature in¬ 
tensity plots [23], we estimate that the splittings at Y is 
about 30 meV at T ~ 20 K. As shown in Fig. [ljf), we 
notice that there is a large hybridization gap between the 
d yz and d xy bands near Y but little hybridization or none 


between the d xy and d xz bands. 

The band structure at Z (Figs. [TJd),[TJg) andQj)) is 
very similar, except for a relative shift along the energy 
direction. In particular, a splitting of about 30 meV is 
observed at Z between the d xz and d yz bands, and an hy¬ 
bridization gap is found between the d xy and d yz bands, 
but not between the d xy and d xz . In Figs. [lJe),[TJ/i) and 
0 *). we show the band structure at M. We distinguish 
two hole-like bands associated with the d yz bands from 
different twin domains. Because of a lack of coherence, 
the d xy electron and hole bands at M are not observed. 
Our data indicate that the splitting at M is about 50 meV 
at T ^ 50 K, which is quite different from the prediction 
of onsite interactions. 

To fully understand the splittings and check if they 
are related, we performed temperature-dependent exper¬ 
iments. The temperature evolution of the d xz /d yz split¬ 
tings at high-symmetry points is illustrated in Fig. [2j 
Except for thermal broadening, the intensity plots show 
that the band dispersions around Y barely change with 
temperature and that the separation between the d xz and 
d yz bands is nearly temperature independent. In other 
words, the d xz /d yz splitting at Y is almost not changed 
within the temperature range studied, and the hybridiza¬ 
tion gap between the d xy band and the d yz band persists 
at high temperature, whereas no hybridization is found 
between the d xy band and the d xz band, indicating that 
none of these phenomena is directly related to the struc¬ 
tural transition. Our conclusion on the splitting at Y 
is reinforced by the comparison of the EDCs at the Y 
point, displayed in Fig. [2] ( u ), and at the Z point (see 
Supp. FigS2). 

Unlike our observation at Y / Z, the band splitting at M 
varies strongly with temperature. The two sets of bands 
from different domains gradually merge with increasing 
temperature. At T = 120 K, we only see one set of 
band structure, which implies the disappearance of the 
domain structure, in agreement with previous results PU¬ 
CE]. The evolution of the EDCs with temperature at M 
is shown in Fig. [2] (v). The dashed lines mark the two 
sets of band tops/bottoms merging at T = 120 K. 

Fig- 0 a) compares the temperature dependence of the 
different splittings and Fig. [3 [b) gives a schematic rep¬ 
resentation of the experimental splittings and hybridiza¬ 
tions observed below and above T s . The splittings at Y 
and Z have the same amplitude, which varies very slowly 
with temperature, even across the structural transition. 
In sharp contrast, the splitting at the M point is nearly 
twice that at the Y point at low temperature, but it de¬ 
creases with temperature and vanishes at 100 - 120 K. 
We conclude that we must introduce two parameters to 
explain the data. The splitting at T/Z is temperature in¬ 
dependent and affects only the BZ area around Y and Z, 
while the parameter inducing the splitting at M only af¬ 
fects the M point and is related to structural transition. 
Since it does not affect the splitting at the BZ center, 
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FIG. 2. (a - e) ARPES intensity plots of the band structure at T at different temperatures. The intensity in each plot is the 
sum of data acquired with C + and C _ polarized photons. (/ - j) EDC curvatures of (a - e). (k - o ) ARPES intensity plots of 
the band structure at the M point at different temperatures, recorded in the a geometry, (p - t) MDC curvatures of (k - o ). 
( u ) EDCs of (a - e) at k x = 0. The blue dashed line is the second derivative of the blue solid line with an extra minus sign. 
The two dashed black lines indicate the peak positions. (T) EDCs of (k - o) at k x = 0. The black dashed lines correspond to 
the EDC peaks. In all the intensity plots, the red lines represent the d yz orbitals, while the blue/cyan ones represent the d xz 
orbitals. All cuts are along T-M or the Z-A high-symmetry direction. All the intensities are divided by Fermi function at the 
corresponding temperatures. 


the order parameter responsible for the splitting at the 
M point must have an anisotropic form of orbital order, 
such as the d-wave orbital order defined on the Fe-Fe 
bonds [20] : 

H-bond — E Am (T) (cos k x cos k y ) (n xz (k) + n yz (k)), 

k 

In Fig §d), we provide detailed calculations and show 
that the d- wave orbital order can explain the experimen¬ 
tal band structure near the M point very well with an 
estimated coupling constant Ao ~ 60 meV in the low- 
temperature limit. 

Two major candidates for the splitting at T/Z are the 
SOC [18j and the onsite ferro-orbital fluctuations |20j. 
However, both explanations contain severe flaws. Indeed, 


SOC can break the glide symmetry that prevents the d xy 
band at k + Q to hybridize with the d xz /d yz bands at k in 
the 1-Fe unit cell (Supp. Part I) |24j[25]. However, such 
hybridization has an equal strength for both d xy ^/d xz ^ 
and d xy ^j d yz ^ hybridizations. Thus, the observation of 
hybridization between the d yz and d xy bands but not be¬ 
tween the d yz and d xy bands is strongly against the SOC 
origin m- In addition, similar splitting at T has been re¬ 
ported to be strongly doping dependent in LiFeAs, which 
is apparently in contradiction with the SOC scenario 
m- The other candidate, the onsite ferro-orbital or¬ 
dering or fluctuations, should remove the d xz /d yz degen¬ 
eracy across the entire momentum space, as illustrated 
by our calculations shown in Figj3|c), which is incon¬ 
sistent with the absence of splitting at M above 120 K. 
Together with the doping-dependent splitting observed 
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FIG. 3. (a) Summary of the d xz /d yz splittings at T, Z and M 

as a function of temperature. The splitting at M disappears at 
a slightly higher temperature than T s , which might be caused 
by short-range ordering or fluctuations above the transition. 
(b) Draft of band structure extracted from experimental data 
and fitted with a tight binding model. The d xy , d xz and d yz 
symbols are for T > T s only since the orbital characters will 
change along T-Mi and T-M 2 . (c - d) Band structure of onsite 
and bond orbital order, calculated from a tight binding model 
(see Supp. Part I for the details). 


in LiFeAs m , we have strong reasons to believe that in 
FeSe the splitting at T and the hybridization between the 
d xy and d yz bands originate from magnetic fluctuations. 
The magnetism in FeSe is more frustrated than in the 
iron-pnictides, and long-range magnetic ordering is thus 
unstable |28l[29]. However, nematicity and magnetic fluc¬ 
tuations can still be strongly coupled [30j[3T]. Thus, the 
splitting at T and the hybridization between d xy and d yz 
observed above T s are very likely signatures of this cou¬ 
pling. In any cases, our current results with two distinct 
dxz/dyz splittings suggest a more complicated interplay 
between the magnetic and orbital degrees in FeSe than 
previously expected. 

In conclusion, we report the temperature evolution of 
the detailed electronic band structure in FeSe single¬ 
crystals. We observe two distinct d xz /d yz band split¬ 
tings at the high-symmetry points. The splitting at M 
is related to the structural transition and has a d -wave 
form factor, while the splitting at T originates most likely 
from magnetic frustration. Our results clearly exclude 
the commonly-believed ferro-orbital order and require a 
new consideration of the origin and implication of the 
orbital order in FeSCs. 
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